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Based on the concepts of artificially microstructured materials, i.e. metamaterials, we present here the 
first practical realization of a radial wave crystal. This type of device was introduced as a theoretical 
proposal in the field of acoustics, and can be briefly defined as a structured medium with radial 
symmetry, where the constitutive parameters are invariant under radial geometrical translations. Our 
practical demonstration is realized in the electromagnetic microwave spectrum, because of the 
equivalence between the wave problems in both fields. A device has been designed, fabricated and 
experimentally characterized. It is able to perform beam shaping of punctual wave sources, and also to 
sense position and frequency of external radiators. Owing to the flexibility offered by the design 
concept, other possible applications are discussed. 

Manipulation of light and its propagation characteristics has been a recurrent problem in the scientific 
community. Recent studies have made use of the so-called metamaterials to modify the propagation 
characteristics of electromagnetic waves almost at will, at least from a theoretical point of view 1,2 . 
These are artificially made microstructures, based on unitary elements of size much smaller than the 
wavelength of the illuminating radiation. The possibility of designing materials' with precise and controllable 
constituent parameters, those that define the solution of the wave propagation equations, has opened a very 
wide range of possibilities going from negative refraction 3 , to sub -wavelength focusing 4 or cloaking 5 . A 
number of proposals exploit the exotic propagation characteristics that can be achieved by the proper design 
of microstructured propagation media 6 " 11 . These concepts can also be extended, with few limitations, to 
different application fields if they share a fundamental wave nature 12 . In this context, the research area has 
become truly multidisciplinary and one can cite as an example the recently proposed Radial Wave Crystals 
(RWCs) 1314 . These microstructured devices are based on the idea of achieving cylindrical structures designed 
in such a way that their constitutive parameters can verify the Bloch theorem. When this is verified, 
propagation of electromagnetic or acoustic waves can be described by means of a dispersion equation or 
a dispersion diagram due to the intrinsic radial periodicity of the microstructures. The fundamental goal of 
these structures is to have a precise control of the allowed and forbidden propagation modes at the different 
frequency bands, since the concept does not present any limitation of frequency scalability. In turn, this 
permits thinking about different types of applications such as resonating elements, beam-shaping shells for 
radiation sources, frequency or position detectors for electromagnetic wave signals or even elements to allow 
a wireless power transfer between them 15 . 

In this article, we present the first practical realization of a Radial Photonic Crystal (RPC), the electromagnetic 
version of the RWC. A 2D device is designed starting from the constitutive parameter functions that determine a 
certain dispersion band profile. This profile is verified in a finite size microstructure implementing a reduced 
version of the ideal configuration, which yields an almost identical behaviour to the full version microstructure. 
The Radial Photonic Crystal is simulated, fabricated and characterized in a measurement setup. Comparison of 
the obtained results is given and analyzed, specially taking into account potential applications. Also, from 
measured and simulated results a parameter extraction is performed in order to recover analytically the order 
of the resonant mode at the frequency of operation. To conclude, we consider various potential applications 
beyond the ones presented in this paper and explain the fundamental aspects of the simulation, fabrication and 
measurement procedures in the Methods section. 
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Results 

Design of a radial photonic crystal. The Radial Photonic Crystal 
design is based on recent theoretical work 13 where it is demonstrated 
that multilayered structures with radially periodic anisotropic 
parameters yield propagation media that verify Bloch's theorem. 
Therefore, it is possible to describe wave propagation by means of 
a dispersion diagram, where the band structure can be tailored 
according to the required objectives. For example, it is possible to 
fix the propagation band of one of the modes of the structure 
completely within the forbidden band of another one, or to fix the 
cut-off frequency of a certain mode coincident with the beginning of 
the forbidden band of a lower order mode. To achieve these and other 
features, constitutive parameters have a radial dependence and are 
defined with tensors: they are anisotropic. We may chose working 
with TM waves to make a practical design, although equivalent 
designs and behaviours exist for TE waves and also for acoustic 
waves. If we take TM Z waves (z-polarized E-field) we have three 
parameters to define: fi r (r), \i Q (r) and s z (r). In view of making a 
practical device we have designed a multilayered structure 
(concentric pairs of alternating layers of types a and b). The 
microstructure has unitary cell parameters, respectively for layers a 
and b, given by the following equations: 
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where d = d a + = 5 + 5 mm is the radial periodicity (see Fig. SI 
for a graphical representation). These functions generate a set of 
three anisotropic parameters that can be implemented by means of 
structured arrays of elements fabricated with standard patterning 
and etching techniques. This system may be scalable in frequency 
by varying periodicity parameter d, with only the obvious limitations 
related to any practical implementation. A radial multilayered 
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Figure 1 | Dispersion diagram of the ideal RPC structure with full set of 
constitutive parameters. Cut-off frequencies of modes with q > 0 are 
equally spaced (aprox. 1.1 GHz). Lower right inset represents a 4 layers 
(finite size) device with internal cavity of radius r = 1 5 mm and periodicity 
d a = d h = 5 mm. Dots and associated E-field patterns correspond to the 
resonance frequencies of the selected operation mode (q = 3;/= 3.8 GHz) 
and adjacent ones of the 4-layers RPC shell. 



structure having these constituent parameters has the dispersion 
diagram shown in Fig. 1. The first 5 modes are represented corre- 
sponding to the unit cell analysis (infinite periodic structure). The 
inset at the lower right corner shows the scheme of a 4 layers struc- 
ture (2 a-type layers and 2 fr-type layers). This finite size structure has 
an internal cavity in order to assure that /n r (r) and e z (r) do not reach 
infinite values at the centre of the shell. Also, superimposed on the 
curves are some of the E-field patterns of the Fabry-Perot resonances 
corresponding to the finite size structure. The symmetry order of 
these modes can be easily identified since the number of lobes in the 
E-field patterns is double with respect to the resonance order. Only 
the first mode (q = 0) does not have a cut-off frequency, the rest have 
almost equally spaced cut-off frequencies (1.1 GHz difference). In 
what follows, the selected operation frequency corresponds to the 
q = 3 Fabry-Perot resonance of the four layers device (f= 3.8 GHz, 
corresponding to a vacuum wavelength k 0 ~ 79 mm). Regard- 
ing the impedance and refractive index values, directly obtained 
from the permittivity and permeability, both have values close 
to the background (vacuum), thus assuring a reasonable match 
between the shell and the surrounding medium (see Fig. S2). 

Although it should be possible to implement a device with the full 
set of parameters previously given, we simplified the structure to a 
reduced parameter set for the sake of fabrication time. The dominant 
parameter in terms of functionality is the angular permeability fie(r), 
therefore it is important to keep it unchanged. It is equivalent to p r 
(the radial component of the inertial mass) in the acoustic counter- 
part of the device 13 . The other two parameters (fi r (r) and s z (r)) are 
replaced by their corrected mean values (see supplementary materi- 
als). With these assumptions, a simple way to implement the per- 
meability and permittivity functions is using arrays of Split Ring 
Resonators (SRRs) 16 , with one metallic ring on each side of a dielec- 
tric substrate. Let us recall that SRRs arrays have a frequency dis- 
persive Lorentz-type resonant magnetic response 17 . This response is 
generated by a magnetic field perpendicular to the rings plane. 
Depending on the geometry of the split rings, it is possible to obtain 
a wide range of potential permeability values at a specific frequency 
(since it is a resonant response giving values from negative to pos- 
itive). A scheme of the proposed structure is displayed in Fig. 2(a). 
This structured resonator array has an anisotropic permeability ten- 
sor with mainly an angular magnetic response in the cylindrical 
coordinate system. Also, it has an electric response derived from 
the dielectric substrate between the metallic rings and the currents 
flowing through the rings. A similar type of structure has been 
employed specifically in a context of beam-shaping applications 18 . 
Our design is formed of 4 concentric arrays (layers) of SRRs having 
the unit cell parameters shown in Fig. 2(b). Exact values of each 
dimension are given in the supplementary materials. The geometric 
parameters that differentiate the concentric layers are the metallic 
pattern widths (w) and the ring gap (g). Each layer is therefore 
formed by an integer number of identical unit cells that synthesise 
the fi 0 (r) values displayed in Fig. 2(c). The variable function of fi 0 (r) 
has a stair-like profile. The continuous variation within each layer 
given by Eq. (2) is here replaced by stepped values of \iq(t) since the 
periodicity parameter d is small enough compared to the operation 
wavelength to permit this discrete approximation. Let us note from 
Fig. 2(c) that \ie(r) values are between 0 and 1, fi r (r) = 1 (no radial 
magnetic response) and s z (r) = 3.4 (constant with the radial coord- 
inate). It is assumed that the fact that the SRRs are organized in 
concentric circular layers does not modify seriously the response 
obtained from the one of a square unit cell as the one used in 
Fig. 2(b). In this case, this approximation is better for the large radius 
or external layers. 

One sample of the fabricated RPC shells is displayed in Fig. 3(a). It 
has an inner cavity with radius r int = 15 mm and an external radius 
°f r ext = 35 mm. The diameter of the shell is smaller than one 
wavelength at the operation frequency. It includes a foam support 
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Figure 2 | (a) Implementation scheme of a 4 layers RPC; (b) unit cell characteristic dimensions and simulation domain for constitutive parameters 
extraction; (c) angular permeability profile as a function of the radial coordinate showing a stepped profile for the 4 layers device. 



of Rohacell (e r ~ 1) to position and hold the SRR arrays. This support 
is neutral to the EM fields. The SRRs are made of metallic strips on 
both sides of dielectric substrate plates. The employed measurement 
setup (see Fig. 3(b)) is a parallel plate waveguide where the RPC is 
illuminated with punctual sources (within the inner cavity or extern- 
ally) or alternatively with a plane wave. Field maps can be obtained 
with a Vector Network Analyzer by varying the position of a sensing 
probe and by recording the transmission and reflection coefficients. 
This is an automated procedure as explained in the Methods section 
that yields magnitude and phase information of the E-fields. 

Characterization. In a first stage, the designed device has been first 
analyzed with a 2D finite element solver 19 by defining the shell 
directly using Equations (1) to (3). Then, constitutive parameters 
have been adapted (and reduced) for the different layers forming 
the device using the values given in Fig. 2(c). In both cases, we call 
this an analytical definition of the RPC. This permits the analysis of 
the multilayered structure in terms of an eigen-value problem to 
extract the Fabry-Perot resonance frequencies displayed in Fig. 1. 
Also, by including point (or line if considering the z-invariance) 
sources, it is possible to obtain field maps of the interaction source 
+ shell. In a second stage, the geometry of the shell can be translated 
to the SRR arrays design. This is performed with a 3D finite element 
solver 20 that allows the optimization of the SRR dimensions to 
produce the target parameters of Fig. 2(c) through standard 



retrieval techniques 21 . These are independent unit cell simulations 
(as from Fig. 2(b)) and the complete procedure for the four layers 
gives a full geometric definition of the RPC. Then, once the geometry 
of each layer is finalized, it is possible to perform numerical 
simulations of the four layer SRR shell (Fig. 2(a)) to obtain the 
field patterns at different frequencies and source positions. These 
are ab-initio simulations that yield the E-field patterns on the top 
plane of the simulation domain of Fig. 2(a), where also the sensing 
probe of the measurement setup captures the measured values. 
Finally, an extensive measurement process is performed for 
different configurations with two prototype devices. Therefore, we 
obtain three sets of comparable results: analytical simulations, 
ab-initio simulations and measurements. 

Analysis of the first configuration consisted of a point source 
located within the inner cavity of the shell (r = 11 mm from its 
centre). This source radiates at a frequency corresponding to the q 
= 3 Fabry-Perot shell resonance (f= 3.8 GHz). Thus, Fig. 4 displays 
the E-field patterns of the combined shell + source structure, where 
the RPC acts as a beam-forming shell. Isotropic radiation of the point 
source is transformed onto a sextupole radiation pattern, corres- 
ponding to the mode order of the excited shell resonance. 
Obviously, this effect can apply to the rest of the resonant modes 
of the shell, at each frequency, with its corresponding beam shaped 
radiation pattern. To excite modes with symmetry q > 0, the inner 
source must be placed at a position displaced from the exact centre of 
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Figure 3 | (a) Radial Photonic Crystal shell including a Rohacell foam 
support for positioning and holding the SRR arrays; (b) measurement 
setup (top plate is removed to allow inner vision) including positioning 
system in two perpendicular axis. 



the cavity. Similar field patterns can be observed on all three panels, 
with slight differences basically due to the discrete implementation 
scheme (this can be seen in panels b) for the ab-initio simulation and 
c) for the experiment). Furthermore, the experiment also shows that 
shell lobes close to the inner source have higher strength. Such an 
effect is more pronounced than in the ideal Comsol simulation. 
These results confirm that, on the one hand, the design procedure 
is validated through the SRR implementation. On the other hand, the 
functionality of the resonant shell is verified, since it can shape the 
radiation pattern of the point source. 

A second alternative is to excite a RPC with an external point 
source at the resonance frequency of the shell. This means that 
depending on the frequency of the external source that illuminates 
the shell, a specific mode can be excited. An example of this is pro- 
vided in Fig. 5, where results from analytical simulation, ab-initio 
simulation and measurement are displayed. The point source is 
located at a distance of 95 mm (> k 0 ) from the centre of the shell 
and it radiates again at the nominal frequency / = 3.8 GHz. If the 
microstructure is designed to resonate at specific frequencies, it can 
detect the presence of radiators at these frequencies that would excite 
the corresponding shell resonances. The shell resonant pattern can 
therefore identify the radiation frequency of the external source 
(within a limited number of frequencies, see Video SI in the supple- 
mentary materials). One additional interesting point of this scenario 
is the fact that the resonant pattern in the RPC shell has one of its 
lobes oriented towards the external source. In this example, it is the 



one with the lowest intensity. This enables locating the position of the 
emitting source according to the resonant pattern orientation. This 
pattern can be especially accurate if several RPCs are used simulta- 
neously at a g = 1 symmetry resonance (see Video S2 in the supple- 
mentary materials). 

For example, if two identical shells are illuminated by the same 
source, then it is possible to use the pair of shells and their respective 
resonance patterns to identify the position of the emitting source. 
Such configuration is depicted in Fig. 5 where two identical shells are 
illuminated by a single point source located approximately at 90 mm 
from their centres. The real part of the E-field is displayed and 
consistent results are observed between all three field maps. The 
main application of this feature is for position detection of point 
sources by triangulation (see the supplementary materials for more 
details on the position determination). At last, a final measurement 
has been performed in order to test the utility of the proposed design 
in view of detecting plane waves (distant sources). In this case, a 
modified setup has been employed with a large parabolic type 
reflector located at the extreme left side of the measurement setup. 
Although imperfect, the reflector produces a plane wave to illumin- 
ate the RPC. Results are displayed in Fig. 7, where it can be noticed 
that the quasi-plane wave again excites the q = 3 resonant mode of 
the shell. The illuminating wave is slightly bended because the dis- 
tance from the reflector to the shell is short and approximately equal 
to two wavelengths. 

As a general comment to these comparisons between simulations 
and measurements, it is easy to observe how the response of first the 
ab-initio simulation (SRR design) and then the fabricated shell are 
progressively degraded with respect to the ideal 2D analytical simu- 
lation. Since ab-initio simulations and measurements are performed 
on a true 3D configuration, only the top plane of the microstructure 
can be mapped. This top plane includes the splits of the rings that 
locally concentrate high E-fields. Although for the 3D cases the 
obtained results clearly reproduce the expected 2D behaviour, this 
effect is superimposed to the shell resonance. 

Extraction of resonant mode coefficients. One additional step in 
the analysis of the measured and simulated results is performed 
through the extraction of resonant mode coefficients from the 
measured or simulated field maps. These field maps represent a 
linear combination of the waves emitted by the source and the 
scattered response of the shell when these waves impinge on it. 
These responses can be evaluated, for cylindrical waves, from the 
summation of Bessel and Hankel functions of different orders (q 
values), with amplitude coefficients that depend on the strength of 
the respective contributing modes. This procedure allows evaluating 
what modes are present, or even dominant, at each frequency and it 
helps to define where the resonance frequencies are located. Starting 
from the configuration of Fig. 5, the measured area can be divided in 
two regions. First region (I) is the one outside of the shell, where the 
direct source contribution and scattered contribution from the shell 
are added. Second region (II) is the inner cavity of the shell where the 
field is only composed of the scattered response from the shell. 
Although this type of analysis can be carried out in both regions 
with the reported difference, we can express the field for the 
second region as: 

*n=Er-coV,(fe-)^. ( 4 ) 

where J q are Bessel functions whose coefficients B q are unknown. A 
linear regression fit is applied to obtain the B q values from the field 
maps of the simulated and measured results. Results of this 
regression analysis are plotted in Fig. 8 for the modes q = 1, 2, 3. 
A high coefficient of determination (R 2 = 0.9986) is obtained from 
the regression analysis of the measured results, which is better than 
for the analysis of the region outside the shell (region I). It is clearly 
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Figure 4 | Real part of the E-field map of a point source placed inside the RPC shell (in the inner cavity) and radiating at a frequency f = 3.8 GHz (q = 3 
mode): (a) analytical simulation; (b) ab-initio simulation; (c) measurement results. Field scale is normalized with respect to the extreme values in each 
plot. 



observed that in a wide frequency range close to the nominal 
resonant frequency the dominant mode has a symmetry order q = 3. 

Discussion 

Several points merit a specific discussion in the design of the pro- 
posed devices. The procedure followed with these microstructures is 
not based on transformation optics, as it is done usually for other 
complex metamaterial- inspired devices 22 . In our case, it is a design 
directly based on dispersion diagram engineering, with selection of 
the allowed and forbidden resonant modes. Essentially the structures 
behave as Fabry-Perot resonators, although inner cavity (Bragg) RPC 
resonances may also be used for creating very high quality factor (Q- 
factor) resonances, as it was demonstrated in acoustics 14 . The par- 
ticular design selected in this case has, on the contrary, very low Q- 
factor resonances. This is mainly due to the fact that the devices 
are specifically adapted to the background. Therefore, impedance 
matching, achieved with reduced impedance values close to one, is 
a target for the microstructures. This fact has a twofold implication. 
On the one hand, it provides quite wide resonances (Q- values are 
close to 10 at the operation frequency). On the other hand, and since 
the resonance is not very sharp, the designed microstructure can 
work not only at the exact resonance frequency, but also over a wide 
range around this nominal value. This means that fabrication and 
measurement uncertainties can be accepted without excessive 



degradation the operation features of the devices. From these char- 
acteristics, the proposed design is well suited for the applications 
described in this work, i.e. interaction with radiation sources. EM 
fields will interact intensely with the shells at and close to the res- 
onance frequencies. 

Another point is related to the slight but systematic difference we 
found between the resonance frequencies of the designed devices, 
and the measured ones. This refers to the fact that simulations in 
Fig. 3 to Fig. 7 are given at 3.8 GHz, whereas measured results are 
actually given at 4 GHz. From the discussion in the previous para- 
graph, we experimentally noted that, good results were obtained at 
3.8 GHz (in terms of agreement with simulations), but the shells had 
better defined resonant patterns at 4 GHz. Although a number of 
factors can have an effect on this problem, we identified the one 
responsible for most of this deviation. Some of the factors can be 
attributed to fabrication tolerances in the geometrical dimensions, 
uncertainties in the exact value of the dielectric constant of the sub- 
strates employed for the SRRs, or imperfections in the design of the 
measurement setup. Systematically we have obtained slightly better 
results at frequencies slightly above the nominal 3.8 GHz frequency. 
This is in a large part due to the fact that our experimental setup is 
based on mobile parallel plates with the RPC in between them. With 
this scheme, it is not possible to exactly separate the two metallic 
plates a distance equal to the RPC shell height (h t ). Otherwise, the 
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Figure 5 | E-field complex magnitude map for a point source illuminating the RPC shell (source - center separation r = 95 mm) and radiating at f = 
3.8 GHz (q = 3 mode): (a) analytical simulation; (b) ab-initio simulation; (c) measurement results. Field scale is normalized with respect to the 
maximum value in each plot. 
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Figure 6 | Real part of the E-field for a point source illuminating 2 RPC shells at a distance of approximately 90 mm to their centres and radiating at a 
frequency f = 3.8 GHz (q = 3 mode): (a) analytical simulation; (b) ab-initio simulation; (c) measurement results. Field scale is normalized with respect 
to the extreme values in each plot. 



relative movement between the plates could easily displace the shell 
and invalidate the measured field map. Let us note also that the 
probing antenna of the top metallic plate is designed to penetrate a 
minimum possible depth inside the inter-plate space. This trade-off 
guarantees a sufficient measured signal level but also a minimum 
invasion in the measured plane (and in principle no contact with the 
shell). In practice approximately a 1 mm air gap was left between the 
top of the shell and the upper metallic plate. This was the minimum 
safety separation to avoid mechanical displacement of the shell with 
the plate movement. Nevertheless, this air layer on top of the shell 
had an impact on the constitutive parameters theoretically synthe- 
sized with the SRR arrays. Actually, it modified them by basically 
increasing the permeability values and decreasing the permittivity 
values. The combined effect produced a slight increase of the res- 
onant frequency, verified by simulating each unit cell (as in Fig. 2(b)) 
including the upper air gap. This explains why results are better at 
4 GHz instead of 3.8 GHz, the nominal resonance frequency. Also, 
in Fig. 8, the maximum of the B 3 coefficient is obtained at around 
4 GHz, although at 3.8 GHz the mode with q = 3 is still the dom- 
inant one. Simulation results in Fig. 8 do include the 1 mm air layer 
(additional simulation details are given in the Methods section). 

Overall, the design flexibility offered by the proposed concept of 
Radial Photonic Crystal allows many different applications for these 
microstructures. At the opposite side of the considered constitutive 



parameters, we may think of a configuration with very high values of 
permittivity and permeability (as high as it could be achieved from an 
implementation point of view). In this complementary scenario, very 
sharp resonances would be obtained with large E-field values con- 
fined within the resonant shells. This situation could be very well suited 
for wireless power transfer applications. In this case, it is intended to 
transfer power between identical resonators that are placed within a 
moderate separation (wireless transfer). The long evanescent fields 
(not radiated as in the reported configuration) may be used for this 
power transfer, since they would provide electromagnetic energy from 
one RPC acting as an active source (or feeder) and a second one acting 
as a receptor. The feeder would be configured with a point source 
within the inner cavity. For this type of application it is mandatory 
to use resonators with high Q values. This is possible by proper engin- 
eering of the constituent parameters so that impedance and refractive 
index allow the devices to highly confine the EM fields within the shell 
region. Energy transfer is then possible through the evanescent mode 
tails outside the shell of the resonance patterns by proximity coupling 
(strong coupling regime). However it is not always necessary that both 
coupled resonators be placed at a very short separation between them 
to make the wireless power transfer efficient 15 . 

In this work, we have demonstrated the feasibility of Radial 
Photonic Crystals, a class of metamaterial-based microstructures 
with a large number of potential applications. We experimentally 




Figure 7 | E-field complex magnitude map for a plane wave illuminating an RPC shell at a frequency f = 3.8 GHz (q = 3 mode): (a) analytical 
simulation; (b) measured results. Results in (a) include a deformed (slightly bended) wave shape according to the practical setup characteristics. Field 
scale is normalized with respect to the maximum value in each plot. 



SCIENTIFIC REPORTS | 2 : 558 | DOI: 1 0.1 038/srep00558 



6 




validated the theoretical concept through the design of a microwave 
shell that permits modifying the radiation characteristics of point 
sources and generating radiation diagrams based on the shell reso- 
nances. Also, these devices permit, by means of a proper design, 
identifying the radiation frequency and position of external point 
sources. An overall good agreement has been obtained between the 
designed shells and the measured results by comparing electric field 
maps in terms of magnitude and real parts. Moreover, the experi- 
mental data also confirm the resonance behaviour by means of an 
analytical extraction of the resonant mode coefficients, and clearly 
identify the resonating mode in the RPC shell. A broad spectrum of 
potential applications is now open for the use of these microstruc- 
tures that are scalable with frequency and can be translated into the 
acoustic domain. 

Methods 

Device fabrication was performed either by standard techniques involving laser 
patterning of a spin- coated polymer layer and chemical etch of the exposed metal in a 
first run or a multistep process combining laser micromachining, protective layer 
definition and chemical etch in a second run. Although the second process guarantees 
an accurate and better definition of all geometrical dimensions, it turns out that the 
first process resulted in a globally better response for the fabricated RPC shell. We 
believe a certain degree of contamination of the composite substrate was generated 
during laser micromachining in the second process, thus degrading the resonant 
response. The substrate used for both processes was Neltec NY9220 with low per- 
mittivity and loss (e r = 2.2 and tan 5 = 0.0009), with uniform dielectric (0.381 mm) 
and metal (35 um) thicknesses for all layers. This observation is in accordance with 
the first point of the Discussion section. Even if there were slight differences in the 
dimensions generated because of the fabrication tolerances of the first process, the 
obtained results in terms of device operation were not affected in a sensible manner. 
On the contrary and as a result of the second fabrication process, metallic particles 
embedded in and contaminating the dielectric substrate did affect more the operation 
of the device. 

The measurement setup consisted of a parallel metallic plate setup, with a mobile 
lower plate with respect to the fixed upper one. The movement between them is linear 
in two orthogonal axes. General details about this type of setup can be found in 
previous works 5 ' 23 ' 24 . Specific characteristics of our setup are described next. The total 
useful measurement area was approximately 280X280 mm 2 , divided in 4 consecu- 
tively measured quadrants. The resolution employed was identical in both movement 
axes, 2 mm, since the sizes of the employed SMA connectors used as emitting and 
sensing probes were about 1 mm in diameter. We used a nominal separation between 



both metallic plates of approximately 10 mm, taking into account that the shell height 
is h t = 9 mm and approximately 1 mm gap is left above the shell to avoid contact with 
the upper plate sensing connector. In the case of the plane wave excitation, a con- 
nector combined with a 2D reflector of a parabolic type is employed as the radiation 
source. The setup allows retrieving information on relative magnitude and phase of 
the electric field. 

The fact of using SRR arrays to implement the RPC shell implies that the dispersive 
behaviour of the SRR unit elements will interact with the expected behaviour of the 
shell as a whole, since constitutive parameters will not be constant with frequency 17 . 
This problem was alleviated by using constitutive parameters in a region where 
dispersion is limited for permeability and almost inexistent for permittivity (see Fig. 
S3 in the supplementary materials). Actually, the synthesized values of \l q are between 
0 and 1, in a region of the Lorentz-type permeability curve where variation with 
frequency is very smooth. This allows limiting the impact of SRR dispersion on the 
operation of the RPC shell. Nevertheless, below 3.5 GHz it is expected that the SRR 
arrays generate negative values of permeability, while approaching the resonance 
frequency, therefore results below that frequency point completely disagree with the 
expected analytical (non-dispersive) RPC behaviour. On the contrary, above the 
operation frequency (3.8 GHz), and since the permeability variation is very smooth, 
the SRR dispersion influence in the RPC shell behaviour is very small and consistent 
results can be obtained from the fabricated shells, as showed in Fig. 8. 

A linear regression analysis is used to obtain the extraction coefficients displayed in 
Fig. 8. This regression analysis is based on identifying the B q coefficients in Equation 
(4) by fitting the measured and simulated field maps to this expression. This permits 
calculating a coefficient of determination R 2 giving the fraction of the field variation 
explained by the equation model. In our case, R 2 is very high for the measurements 
because we use the field map only inside the inner cavity. This limits the external 
sources of error that impact more on the field pattern external to the RPC. It is also 
possible, from the residual variance of the measured data, to calculate the 95% con- 
fidence intervals as plotted in Fig. 8. Confidence intervals for the simulated data are 
not included in this plot, since for the simulations R 2 is equal to 1 in practice. 

Finally, the extracted coefficients in Fig. 8 are directly obtained from the measured 
field map by processing it. In the case of the ones derived from the analytical simu- 
lation, they require some additional explanation. Taking into account the slight 
resonance frequency shift generated by the air layer above the RPC, we performed a 
new series of simulations of the unit cells (as in Fig. 2(b)) but this time including an 
extra 1 mm air layer above the SRR. This generated a set of new constituent para- 
meters at 3.8 GHz (the operation frequency). This variation is at the origin of the 
observed frequency shift. In order to check that this agreed with the experimental 
observation, we performed a new analytical simulation with the second set of con- 
stitutive parameters (those obtained at 3.8 GHz and invariant with frequency) to 
obtain a new field map for the RPC illuminated by the external source, as in the 
measurement. This time the analytical resonance frequency of the q = 3 mode is also 
shifted upwards and in practice matches the one obtained in the measurements, as it is 
confirmed in Fig. 8. 
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